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SECTION  1.0  GENERAL 


1.1  Purpose  of  the  User's  Manual,  This  user's  manual  for  the  wind  statistics  program  (DNDWST) 
provides  the  information  necessary  to  effectively  use  the  program. 

1.2  Project  References. 

1.2.1  Project  Request,  the  Space  and  Missile  Test  and  Evaluation  Center  (SAMTEC)  sent  a  project 
request  to  the  U5AF  Environmental  Technical  Applications  Center  (USAFETAC)  on  10  March  1978 
describing  the  work  required  in  accordance  with  AWSR  105-18.  The  USAFETAC  project  number  assigned 
for  the  development  of  this  program  was  1893,  task  01. 

1.2.2  Documentation  Concerning  Related  Projects.  NASA  TO  X-73319,  ’Vector  Wind  and  Vector  Wind 
Shear  Models,*  by  0.  E.  Smith,  contains  a  description  of  the  equations  used  in  the  program  and 
outlines  the  theory  involved. 

1.2.3  Documentation  Standards  and  Specifications.  The  documentation  complies  with  DOD  Standard 
7435. 1-S,  Automated  Data  Systems  Documentation  Standards,  dated  13  September  1977. 

1.2.4  Programming  Conventions.  American  National  Standards  Institute  (ANSI)  FORTRAN  programming 
conventions  have  been  adhered  to,  except  in  the  following  cases: 

a.  Literals  within  FORMAT  statements  are  written  inside  apostrophes  instead  of  using  the  "H* 
(Hollerith)  format  code. 

b.  The  integer  variable  ’ITERM",  which  is  always  set  to  -1,  is  the  FORTRAN  logical  unit  for  all 
input  and  output.  This  feature  allows  proper  transfer  of  data  across  the  ARPA  network. 

c.  Free  field  formats,  such  as  "21"  for  two  integers  and  "IF"  for  one  real  number,  are 
occasionally  used  to  give  the  user  more  flexibility  when  typing  an  input  line. 

1.3  Terms  and  Abbreviations 


ARPA:  Advanced  Research  Projects  Agency 
BB&N:  Bolt,  Beranek  &  Newman,  Inc. 

Bivariate  Normal:  A  commonly  used  statistical  distribution  which  is  specified  by  the  means, 
standard  deviations,  and  correlation  coefficient  of  two  normally  distributed  variables. 

DEC:  Digital  Equipment  Corporation. 

Meridional  Component:  The  north/south  component  of  the  wind. 

Percentile  Wind  Speed:  A  wind  speed  which  is  greater  than  or  equal  to  the  actual  wind  speed  for  a 
given  percentage  of  observations.  For  example,  if  a  95th  percentile  wind  speed  is  50  knots,  one 
would  expect  95  percent  of  the  observed  wind  speeds  to  be  less  them  or  equal  to  50  knots. 

Scalar  Wind  Speed:  The  observed  wind  speed,  i.e. ,  how  fast  the  air  is  actually  moving. 

Wind  Speed  Units:  Any  type  of  wind  speed  can  be  used  with  this  program.  However,  the  units  used 
for  the  wind  speed  and  the  velocity  must  be  consistent  throughout  the  program;  for  example,  if  the 
standard  deviations  are  in  meters  per  second  (MPS) ,  then  the  u  and  v  components  must  be  in  MPS,  and 
a  request  for  a  wind  speed  range  will  also  use  MPS. 

Zonal  Component:  The  east/west  component  of  the  wind. 

1.4  Security  and  Privacy.  This  manual  and  the  program  that  it  describes  are  unclassified  and  can 
be  released  to  the  public.  No  privacy  restrictions  are  associated  with  the  use  of  this  program. 
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SECTION  2.0  SYSTEM  SIMIARY 


2.1  System  Application,  The  wind  statistics  program  provides  the  user  with  information  relating  to 
the  probabilities  of  various  ranges  of  wind  directions  and  wind  speeds,  using  an  assumed  bivariate 
normal  distribution  of  the  zonal  find  meridional  components  of  the  wind.  It  allows  the  user  to 
obtain  probabilities  of  extreme  event  ranges  (for  example,  extremely  high  or  low  wind  speeds)  which 
cannot  be  computed  easily  and  accurately  by  direct  access  to  meteorological  data.  This  program  can 
also  recompute  a  set  of  basic  wind  statistics  using  a  rotated  set  of  coordinates,  such  as  cartesian 
coordinates  with  one  axis  oriented  along  the  azimuth  of  a  flight  path. 

2.2  System  Operation.  Users  (who  are  typically  USAFETAC  analysts)  must  first  obtain  a  set  of  five 
basic  wind  statistics  (means  and  standard  deviations  of  the  zonal  and  meridional  components  and 
their  correlation  coefficient)  for  the  station (s)  of  interest  before  accessing  the  program.  The 
"Uniform  Summary  of  Winds  Aloft  Observations"  data  summaries,  produced  by  USAFETAC,  contain  these 
basic  statistics  for  several  hundred  upper-air  observation  stations,  calculated  by  month  of  the  year 
and  season  for  various  altitudes.  These  statistics  are  then  input  to  the  program,  after  which  the 
appropriate  choice  of  option (s)  causes  the  program  to  print  the  desired  information.  The  output  may 
also  be  generated  in  a  form  suitable  for  transfer  across  the  ARPA  network  or  on  tape.  The  overall 
functional  diagram  is  shown  schematically  in  Pigure  1, 


Figure  1.  Overall  Functional  Diagram. 


2.3  Program  Configuration.  The  wind  statistics  program  was  designed  for  a  DEC-10  computer.  It 
presently  operates  on  Bolt,  Beranek  and  Newman’s  system  "B"  (B8NB)  on  the  ARPA  network.  Changing 
the  nonstandard  items  listed  in  Section  1.2.4  of  this  manual  to  ANSI  Standard  FORTRAN  will 
effectively  make  this  program  conpatible  with  other  FORTRAN  compilers. 

The  configuration  necessary  for  a  computer  to  use  this  program  typically  requires  a  FORTRAN 
compiler,  linking  loader,  and  interactive  real-time  input/output  capability. 

2.4  Program  Organization.  The  program  consists  of  a  main  program,  seven  subroutines,  and  eight 
functions. The  main  program  solicits  and  controls  all  input  data,  and  checks  this  data  for  obvious 
errors.  Most  of  the  output  to  the  user's  terminal  or  CRT  is  also  controlled  by  the  main  program. 
The  numerical  calculations  are  performed  by  the  various  functions  and  subroutines;  the  number  and 
name  of  subprograms  involved  depends  on  the  type  of  output  requested  by  the  user.  Same  of  the 
subroutines  and  functions  are  accessed  indirectly  by  other  subprograms.  The  calling  structure  of 
the  main  program  and  subprograms  is  illustrated  in  Figure  2.  The  purpose  of  each  subprogram  is 
outlined  in  Table  1. 
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TABLE  1.  PURPOSE  OF  THE  SUBPROGRAMS 


Subprogram  Purpose 

Subroutine  EL IPSE  Calculates  and  prints  sets  of  points  located  on  an  elliptical  probability 

contour . 

Subroutine  PRSRD  Computes  the  probability  of  the  wind  being  within  specified  range  of  speeds 

and  a  specified  range  of  directions. 

Subroutine  ROTATE  Calculates  a  new  set  of  basic  wind  statistics  based  on  a  rotation  of  the 

X-Y  axes  through  a  given  angle. 

Subroutine  RSPGDR  Gives  the  (conditional)  probability  of  a  specified  range  of  wind  speeds 

when  the  wind  direction  is  known. 

Computes  the  (conditional)  probability  that  the  wind  speed  is  less  than  or 
equal  to  a  specified  value  at  a  single  wind  direction. 

Calculates  percentile  wind  speeds  by  initially  estimating  an  answer  and 
then  using  the  method  of  half  intervals  to  get  a  final  result. 

Calculates  and  prints  a  table  of  probabilities  of  various  ranges  of  wind 
directions  and  speeds. 

These  functions  correspond  to  the  terms  "A",  "B",  "C",  and  "D“  in  Equation 
(32),  (33),  and  (36)  of  NASA  TM  X-73319.  Each  of  these  terms  is  coded  as  a 
separate  function  to  insure  compatability  with  this  document. 

Gives  the  probability  density  function  for  a  specific  wind  direction. 

Confutes  the  probability  density  for  the  normal  distribution  at  a  specified 
point. 

Function  PRDIR  Confutes  the  probability  of  a  specified  range  of  wind  directions. 

Function  TRANS  Modifies  the  input  wind  direction (s)  so  that  all  wind  velocity  vectors 

point  away  from  the  origin. 

2.5  Performance 

2.5.1  Overall  Performance  Capabilities.  The  wind  statistics  program  can  give  the  user  the 
probability  of  a  wide  variety  of  wind  speed  and/or  direction  ranges.  These  ranges  are  subject  to 
certain  limitations  imposed  by  the  internal  accuracy  of  the  program,  irput/output  formats,  and  the 
numerical  approximations  used  in  calculating  certain  types  of  probabilities.  These  limitations  are 
discussed  in  greater  detail  in  Section  2.5.3.  All  users  of  this  program  should  also  remember  that 
the  program  uses  an  assumed  bivariate  normal  distribution  of  the  zonal  and  meridional  wind 
components  to  calculate  its  results.  The  accuracy  of  the  data  output  by  this  program  depends  to  a 
large  extent  on  how  well  the  statistics  of  the  actual  wind  "fit"  the  bivariate  normai  distribution. 
The  accuracy  of  this  distribution's  fit  varies  with  altitude  and  location;  this  variation,  and  some 
methods  of  determining  how  accurate  this  distribution  is,  will  be  dealt  with  in  detail  in  a  future 
USAFETAC  technical  note. 

2.5.2  Response  Time.  The  response  time  of  this  program  depends  on  the  processing  time  required  by 
the  option  (s)  being  used  and  also  on  the  number  of  users  concurrently  running  other  programs  on 
BBN's  system-B.  If  the  usage  of  BBN's  system-B  is  unusually  heavy,  all  the  printout  produced  by  the 
program  will  be  slowed  and  delays  up  to  several  minutes  in  length  may  occur.  Waiting  times  are 
generally  proportional  to  the  CPU  time  required  by  an  individual  program  option.  The  response  time 
of  any  given  option  will  be  several  times  longer  during  periods  of  heavy  BBNB  usage.  The  time  lag 
between  user  input  and  program  output  may  approach  5  minutes  for  some  options  during  periods  of 
extremely  heavy  machine  use.  The  CPU  times  required  by  the  different  options  are  discussed  in 
Section  2.5.5. 

2.5.3  Input/CXitput  Limitations.  The  input  and  output  limitations  existing  in  the  wind  statistics 
program  fall  into  two  categories:  limitations  due  to  formats,  and  truncation  errors  inherent  in  the 
calculations  producing  the  output.  The  limitations  which  apply  to  specific  options  are  listed 
below. 


Subroutine  3LBGDR 

Subroutine  SPCNT 

Subroutine  WROSE 

Functions  A,  B,  C,  D 

Function  PDIR 
Function  PNORM 
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a.  The  permissible  range  of  output  data  points  in  option  #1  (points  on  an  elliptical 
probability  contour)  is  -99.99  to  999.99  wind  speed  units.  Data  values  with  absolute  values  less 
than  0.001  are  shown  as  zero.  Identical  restrictions  apply  to  the  extreme  values  of  x  and  y  printed 
by  this  option. 

b.  All  probabilities  output  by  options  2,  3,  4,  and  5  are  limited  to  the  range  0<x<l  i Aere  x  is 
the  output  probability  value.  The  probabilities  whose  absolute  value  is  less  them  0.0001  are  shown 
as  zero. 

c.  The  wind  speeds  used  as  input  to  options  3,  4,  and  5  should  be  within  the  range  of  0  to  500 
wind  speed  units.  Values  outside  this  range  may  cause  unpredictable  error  conditions  and  abort  the 
program. 

d.  Hie  basic  wind  statistics  output  by  option  #6  (involving  rotation  of  coordinates)  have 
several  limitations.  Hie  mean  wind  components  printed  must  be  smaller  than  10^  units  and  larger 
than  10~4  units;  wind  component  values  with  an  absolute  value  of  less  than  0.01  units  are  printed  as 
zero.  The  correlation  coefficient  printed  appears  as  zero  if  its  absolute  value  is  less  than  0.001. 

e.  Wind  speeds  output  by  option  #8  must  be  smaller  than  10^.  Speeds  smaller  than  0.001  are 
printed  as  zero. 

f.  Probabilities  outfxit  by  option  #9  are  affected  by  the  restriction  given  in  a.  above  with  one 
exception:  the  number  output  for  the  total  probability  may  be  slightly  larger  than  1.  This  number 
will  almost  always  be  between  0.999  and  1.001. 

Several  general  limitations  exist  in  the  program.  These  limitations  apply  in  addition  to  the 
limitations  given  for  specific  options.  The  probability  of  a  range  of  directions  is  calculated 
internally  in  the  program  to  a  precision  of  10“^,  while  the  probability  of  a  range  of  speeds  or  a 
range  of  speeds  and  directions  is  calculated  to  a  precision  of  10-4. 

2.5.4  Error  Detection  Capabilities.  The  wind  statistics  program  is  capable  of  detecting  several 
types  of  input  errors.  The  program  checks  all  input  speeds  for  negative  values;  values  less  than 
zero  produce  the  message 

"ONE  OR  MORE  OF  THE  INPOT  SPEEDS  ARE  NEGATIVE" 

after  which  the  user  is  asked  to  input  the  speeds  again.  Similar  consequences  result  when  a  user 
tries  to  input  a  negative  standard  deviation  or  a  correlation  coefficient  which  is  outside  its  legal 
range  (-1  to  +1).  An  input  percentile  value  for  option  #8  which  is  negative,  zero,  or  larger  than 
or  equal  to  100  will  also  produce  a  warning  message  and  a  request  to  input  the  data  again.  The 
program  will  reply  to  an  illegal  option  number  (less  than  0  or  greater  than  9  by  asking  the  user  if 
the  program  should  be  stopped.  A  "yes"  reply  stops  the  program;  a  reply  of  "no"  allows  the  user  to 
input  a  new  option  number. 

2.5.5  Processing  Times.  The  CPU  time  used  by  this  program  depends  primarily  on  the  output  options 
used  by  the  analyst.  Hie  total  CPU  time  for  a  given  run  can  be  calculated  aproximately  by  adding 
the  time  required  to  print  error  messages,  the  time  required  to  start  and  stop  the  program,  and  the 
time  required  for  each  option  used.  All  CPU  times  given  are  valid  for  a  DEC-10  computer. 

The  CPU  time  required  to  start  the  program,  input  a  set  of  basic  statistics,  and  stop  it  is  about 
0.5  sec.  This  represents  unavoidable  "overhead”  time.  An  additional  0.9  sec  is  required  to  print 
the  optional  informative  summary  of  the  various  available  output  options.  Each  error  message  also 
requires  about  0.2  seconds  to  print. 

Typical  CPU  times  for  the  various  options  are  as  follows: 

a.  Cation  #1  (points  on  an  elliptical  probability  contour)  has  CPU  requirements  which  depend  on 
the  number  of  points  generated.  Ten  sets  of  points  require  about  0.8  seconds,  50  sets  require  about 
2.3  seconds,  and  100  sets  require  about  4.5  seconds. 

b.  Option  ♦ 2  (the  probability  of  a  range  of  wind  directions)  requires  about  0.6  seconds  of  CPU 
time. 

c.  Option  #3  (the  probability  of  range  wind  speeds)  requires  anywhere  from  2  to  15  seconds  of 
CPU  time.  Small  ranges  of  speeds  are  processed  relatively  rapidly. 
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d.  Option  #4  (the  probability  of  a  range  of  wind  speeds  and  directions)  uses  from  4  to  IS 
seconds  of  CPU  time,  depending  on  the  size  of  the  range  involved. 

e.  Option  #5  (the  probability  of  a  range  of  wind  speeds  given  a  wind  direction)  uses  about  0.5 
seconds  of  CPU  time. 

f.  Option  *6  (rotation  of  the  coordinate  system)  uses  approximately  0.2  seconds  of  CPU  time. 

g.  Option  #7  (input  of  new  basic  parameters)  requires  0.2  seconds  of  CPU  time. 

h.  Option  #8  (percentile  windspeeds)  uses  larger  amounts  of  CPU  time.  Typical  CPU  times  for 
this  option  ran  ,e  from  12  to  25  seconds. 

i.  Option  #9  (probability  tables  for  wind  direction  and  speed  ranges)  is  the  most  time 
consuming  option  in  this  program.  Typical  CPU  times  for  this  option  range  from  four  to  six  minutes. 

2.5.6  Flexibility.  The  program  is  designed  to  allow  additional  output  options  to  be  included  at  a 
later  date.  For  example,  sane  possible  additional  options  include  synthetic  vector  wind  profiles, 
vector  wind  shear  probabilities  and  the  modulus  of  the  vector  wind  shear.  These  additional  options 
would  use  some  of  the  subprograms  already  coded  in  the  programs  and  "build"  on  their  "foundation." 
Each  output  option  is  presently  controlled  by  one  segment  of  the  main  program;  additional  options 
would  require  additional  segments  of  code  in  the  main  program  and  one  or  more  additional 
subroutines. 

2.6  Data  Base.  None  required  for  this  program. 

2.7  General  Description  of  Inputs,  Processing,  and  Outputs 

2.7.1  Inputs.  All  input  to  the  program  is  typed  on  a  terminal  or  CRT.  Hie  user  must  have  the  five 
basic  wind  statistics  (mean  zonal  and  meridional  wind  components,  their  standard  deviations,  and 
correlation  coefficient)  available  for  input  to  the  program  before  running  it.  The  users,  typically 
USAFETAC  analysts,  should  also  determine  in  advance  which  option (s)  are  appropriate  for  their  work. 
For  example,  option  #3  would  be  appropriate  to  answer  the  question,  "what  is  the  probability  of  a 
wind  speed  greater  them  20  knots?"  A  list  of  output  options,  option  numbers,  and  descriptions  of 
the  resulting  output  are  given  in  Section  3.2.  The  user  should  then  get  the  reference  information 
required  for  these  option (s);  for  example,  the  number  of  sets  of  points  to  be  given  by  option  #1  or 
the  wind  direction  range  for  option  #2.  Details  on  the  type  and  nature  of  reference  information 
required  for  each  option  are  given  in  Sections  3.2.1  and  3.2.3. 

2.7.2  Processing.  The  user  begins  by  typing  a  set  of  five  basic  wind  statistics.  The  program  will 
then  ask  the  user  for  an  option  number.  After  the  user  types  the  appropriate  digit,  the  program 
asks  for  reference  information.  The  calculation  and  printing  of  output  takes  place  after  all  the 
reference  data  has  been  input.  The  program  then  asks  for  another  option  number.  The  user  can  then 
stop  the  program  (option  #0) ,  change  the  basic  wind  statistics  (option  #7) ,  or  use  another  option. 
The  program  continues  in  this  fashion  until  the  user  inputs  an  illegal  option  number  and  indicates 
that  a  program  stop  is  desired,  or  until  the  program  is  stopped  by  option  #0.  This  process  is  shown 
in  Appendix  A. 

2.7.3  Output.  All  output  is  normally  sent  to  the  user's  input  device  (a  terminal  or  CRT) .  The 
format,  type,  and  content  of  the  output  will  vary  depending  on  the  option  selected.  The  types  and 
formats  of  output  are  described  in  detail  in  Section  3.3.1. 
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SECTION  I  STAFF  FUNCTIONS  RELATED  TO  TECHNICAL  OPERATIONS 


3,1  lnitiati.cn  Procedures.  The  user  should  l  >)  into  'he  ETAC-ON  direct'''"  .. .  BBNB  on  the  ARPA 
network.  USAFETAC  users  requiring  assistance  with  ARPA  acre’-'-  t.aoedures  should  refer  to  the 
USAFETAC  ARPA  user's  guide.  Typira  "DNDWPT"  af*-'"  c  prompt  will  start  the  proqram.  The 

program  will  reply  with 


"***  URAFETAC'/DND  WIND  STATISTICS  PROGRAM  ***" 
to  indicate  that  the  program  has  begun  execution. 

3.2  Staff  Input  Requirements.  The  user  inputs  data  interactively  to  this  proqram  via  a  terminal  or 
CRT.  The  program  always  explicitly  asks  for  any  input  required  at  the  appropriate  times  during  its 
execution. 

3.2.1  Input  Contents  and  Format.  'Che  format  of  data  input  to  the  program  depends  on  the  type  of 
data  the  program  is  requesting.  After  the  printout  oi  its  title,  the  program  will  ask  the  user  if 
he  wants  information  on  the  available  output  options.  A  “Y”  input  produces  a  list  of  these  options, 
after  which  the  program  asks  for  a  set  of  five  basic  wind  statistics.  Typing  anything  other  than  a 
"Y"  allows  the  program  to  go  directly  to  this  .step.  The  basic  wind  statistics  must  be  typed  on  one 
line  in  the  following  order: 

a.  The  mean  zonal  component  of  the  wind,  must  be  typed  as  a  real  nunber  (with  a  decimal  point)  . 
Leading  zeros  may  be  omitted. 

b.  The  standard  deviation  of  the  zonal  component  of  the  wind,  "STDEVX,"  must  be  a  positive  real 
number . 

c.  The  mean  meridional  component  of  the  wind,  must  be  a  real  number. 

d.  The  standard  deviation  of  the  meridional  component  of  the  wind,  "STDEVY,"  must  be  a  positive 

real  number. 

e.  The  correlation  coefficient  of  the  two  components  of  the  wind,  indicated  as  "CORR.  COEFF." 
in  the  program's  request  for  input,  must  be  a  real  number  between  -1  and  1. 

These  five  values  must  be  typed  in  separated  by  comnas  and  without  intervening  blanks.  The  total 
number  of  digits  in  any  input  value  should  not  exceed  eight;  any  excess  digits  will  be  truncated 
after  input  to  the  program  is  complete. 

The  program  asks  for  an  option  number  after  the  basic  statistics  have  been  input.  This  number 

should  be  a  single  digit  from  0  to  9.  The  proqram  will  accept  larger  numbers  and  give  a  warning 

message;  the  use  of  a  decimal  nurttoer  will  give  erratic  results.  The  input  contents  and  format  used 
with  the  individual  program  options  are  shorn  in  Table  2. 

TABLE  2.  INPUTS  TO  THE  PROGRAM  OPTIONS. 

Option  Number  Type  of  Input 

1  a.  A  decimal  probability  value  between  0  and  1 
b.  The  nunber  of  points  required  (an  integer) 

2  TWo  decimal  numbers  indicating  a  range  of  directions 

3  Two  decimal  numbers  indicating  a  range  of  speeds 

4  a.  Two  decimal  numbers  indicating  a  range  of  directions 
b.  Two  decimal  numbers  indicating  a  range  of  speeds 

5  a.  A  single  decimal  number  indicating  a  direction 
b.  Two  decimal  numbers  indicating  a  range  of  speeds 

6  A  decimal  number  indicating  an  angle  of  rotation 

7  Five  decimal  numbers  (basic  wind  statistics).  The 
second  and  fourth  numbers  (standard  deviations)  must  be 
positive;  the  fifth  number  (correlation  coefficient) 
must  be  between  -1  and  1. 


8  A  decimal  percentile  value  between  0  and  100. 

9  a.  Four  sets  of  alphanumeric  characters  describing  the 
period  of  record,  layer  or  level,  location,  and  wind 
speed  units  associated  with  the  probability  table 
being  printed.  The  maximum  lengths  of  these  four 
character  strings  are  20,  30,  25,  and  10  characters, 
respectively. 

b.  A  character  (Y  or  N)  indicating  the  users  choice 
of  where  the  output  should  go.  If  the  user  replies 
with  "Y",  two  more  inputs  are  required: 

c.  A  set  of  up  to  five  alphanumeric  characters 
indicating  the  file  name  to  be  used  for  the  output 
produced  by  this  option,  and 

d.  Another  character  (Y  or  N)  indicating  whether  or 

not  the  program  should  stop  after  completing  this  option. 

3.2.2  Composition  Rules.  Hie  programming  language  used  (FORTRAN) ,  together  with  the  irput  formats. 
Impose  certain  limits  on  the  composition  of  input.  These  limits  are  as  follows: 

a.  Numbers  input  to  the  program  may  not  include  embedded  blanks. 

b.  Hie  characters  which  may  be  typed  in  reply  to  a  request  for  numerical  data  are  limited  to 
the  digits  0-9,  periods,  commas,  the  plus  and  minus  signs.  Typing  any  other  character  will  abort 
the  program. 

3.2.3  Input  Vocabulary.  Only  one  set  of  codes  is  used  in  the  wind  statistics  program.  Hie  output 
options  are  designated  by  numbers  from  0  to  9.  These  options  are  listed  in  Table  3. 

TABLE  3.  PROGRAM  OPTICUS 

Option  Number  Description  of  the  Resulting  Output 

0  Stops  the  program 

1  Prints  points  defining  an  elliptical  probability 
contour.  A  wind  vector  with  the  given  basic  wind 
statistics  will  have  a  probability,  equal  to  a 
value  specified  in  the  user's  input,  of  being 
inside  this  elliptical  contour. 

2  Prints  the  probability  of  a  range  of  wind 
directions.  Hie  range  is  always  taken  to  be 
going  clockwise  from  the  first  direction  given 

to  the  second  one.  Directions  are  always  in  degrees. 

3  Prints  the  probability  of  a  range  of  speeds.  Hie 
wind  speed  units  used  must  be  identical  to  those 
used  for  the  five  basic  wind  statistics.  Hie 
probability  of  a  wind  speed  outside  the  input 
range  will  be  given  if  the  larger  of  the  two  speeds 
defining  the  range  is  input  first. 

4  Prints  the  joint  probability  of  a  range  of  directions 
and  range  of  speeds.  The  conventions  used  in 
determining  ranges  are  identical  to  those  listed  for 
options  2  and  3. 

5  Prints  the  conditional  probability  of  a  range  of  wind 
speeds  when  the  wind  direction  is  known. 

6  Rotates  the  coordinate  system  used  for  the  basic  wind 
statistics  and  prints  the  new  basic  statistics. 

These  new  values  are  then  used  in  all  program 
calculations  until  changed  again  by  the  use  of 
options  6  or  7. 
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Allows  the  user  to  input  a  new  set  of  basic  wind 
statistics.  The  new  statistics  are  then  used  in  all 
subsequent  program  processing  until  they  are  changed 
again  by  options  6  or  7. 

Prints  wind  speed  percentiles.  The  wind  speed  will 
have  a  probability,  equal  to  a  value  input  by  the  user, 
of  being  smaller  than  or  equal  to  the  wind  speed  printed 
by  this  option. 

Prints  a  table  containing  the  probabilities  of  various 
wind  speed  and  direction  intervals.  This  table  is 
labeled  with  the  station  location,  level,  period  of 
record,  and  wind  speed  units  input  by  the  user.  The 
wind  direction  ranges  used  by  this  option  correspond 
with,  and  are  centered  on,  the  standard  directions 
in  a  wind  rose:  N,  NNE,  NE,  ENE,  etc.  The  wind 
speeds  are  broken  down  into  ranges  of  less  than  1 
unit  (calm) ,  1-10  units,  10-20,  20-30,  30-40, 

40-50,  50-60,  60-75,  75-100,  100-150,  150-200,  and 
greater  than  200  units. 


3.2,4  Saaple  Inputs.  A  sanple  user  reply  to  a  program  request  for  basic  statistics 
"INPUT  MEAN  X,STDEVX,MEAN  Y,STDEVY,  CORR.  COEFF."  would  be  "1.3, 3. 4, -22. 7, 6. 8, .2312" 


where 


1.3  =  mean  tonal  wind  conponent 

3.4  =  standard  deviation  of  the  zonal  wind  conponent 
-22.7  *?.  the  mean  meridional  wind  component 

6.8  »  standard  deviation  of  the  meridional  wind  component 
.2312  *  the  correlation  coefficient  of  the  two  wind  components 


The  type  and  format  of  other  inputs  to  the  program  depend  on  the  option  being  used.  Some  sanple 
inputs  for  each  option  which  asks  for  input  are  given  below: 

a.  Sample  input  for  option  #1.  The  user's  reply  to  this  option's  initial  request  "INPUT 
PROBABILITY  VALUE"  is  .23  indicating  that  the  23  percent  probability  contour  is  needed.  The  next 
request  is  "INPUT  THE  NUMBER  OF  X-YLOW-YHIGH  SETS  DESIRED."  Typing  "25"  asks  this  option  to  print 
out  25  sets  of  points. 


b.  Sample  input  for  option  #2.  The  user's  reply  to  this  option's  request  "INPUT  RANGE  OF 
DIRECTIONS"  is  "125.3,178.6"  which  specifies  the  range  of  directions  from  125.3°  to  178.6°.  Typing 
"330., 30."  instead  would  specify  the  range  of  directions  from  NfW  through  N  to  NNE. 


c.  Sample  input  for  option  *3.  Typing  "3., 6."  in  response  to  this  option's  request  "INPUT 
RANGE  OF  SPEEDS,"  tells  this  option  to  get  the  probability  of  a  wind  speed  between  3  and  6  wind 
speed  units.  Typing  in  "10., 0."  instead,  asks  this  option  to  generate  the  probability  of  a  wind 
speed  greater  than  10  wind  speed  units  (corresponding  to  the  range  outside  0-10  units.) 

d.  Sample  input  for  option  #4.  This  option  solicits  input  with  the  statement  "INPUT  RANGE  OF 
DIRECTIONS."  The  user  types  in  "30., 63. 5".  This  option  then  asks  for  more  input  by  printing  out 
"INPUT  RANGE  OF  SPEEDS."  The  user  then  types  in  "0.,33.5".  The  program  will  then  calculate  the 
joint  probability  of  a  wind  direction  between  30  and  63.5  degrees  (NNE  to  approx.  ENE)  and  a  wind 
speed  below  33.5  units. 


e.  Sample  input  for  option  #5.  This  option  prints  "INPUT  DIRECTION."  The  user  then  types  in 
"225.".  This  option  then  prints  "INPUT  RANGE  OF  SPEEDS."  The  user  then  types  in  "0.,33.5.".  The 
program  will  then  calculate  the  probability  that  the  wind  speed  will  be  equal  to  or  less  than  33.5 
units  when  the  wind  direction  is  225  degrees  (SW) . 


f.  Sample  input  for  option  16.  This  option  asks  for  input  with  the  statement  "INPUT  ANGIE  OF 
ROTATION."  Typing  "15."  will  indicate  to  the  program  that  a  counterclockwise  rotation  of  15°  is 
required. 


g.  Sanple  input  for  option  #8.  This  option  states,  "INPUT  TOE  REQUIRED  PERCEWTIIE."  Typing 
"88.0"  asks  the  program  to  compute  the  88th  wind  speed  percentile  value. 
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h.  Sample  input  for  option  IS.  This  option  requires  several  lines  of  input.  The  lines 
enclosed  in  quotes,  below,  are  entries  typed  by  the  user.  The  other  lines  are  typed  by  the  program. 

INPUT  TOE  PERIOD  OF  RECORD  -  UP  TO  20  CHARACTERS 
"JAN  70  -  DEC  74" 

INPUT  TOE  IEVEL  OR  LAYER  -  UP  TO  30  CHARACTERS 
"4500  FEET  -  5000  FEET" 

INPOT  TOE  LOCATION  -  UP  TO  25  CHARACTERS 
"ANDREWS  AFB" 

INPOT  THE  WIND  SPEED  UNITS  -  UP  TO  10  CHARACTERS 

"MFH" 

DO  YOU  WM7T  TOE  PROBABILITY  TABLE  WRITTEN  OUT  TO  A  FILE?  (Y*=YES,  N»NO) 

Nyn 

PLEASE  INPOT  TOE  FILENAME  -  UP  TO  5  CHARACTERS 
"NNNNN" 

WOULD  YOU  LIKE  TOE  PROGRAM  TO  STOP  AFTEIWARD?  (Y=YES,  N=NO) 

"N" 

YOUR  PROBABILITY  TABLE  WILL  BE  SENT  TO  FILE  NNNNN.DAT 

Option  #7  asks  for  another  set  of  basic  statistic-,  while  option  #0  does  not  require  input. 

3.3  Output  Requirements.  All  the  output  appears  on  the  user's  input  device  (normally  a  CRT  or 
terminal).  The  output  varies  depending  on  the  option  nuntoer  selected.  Option  #1  generates  a  set  of 
points  for  plotting  a  probability  contour.  Options  2,  3,  4,  and  5  generate  probabilities  of  wind 
speed  and/or  wind  direction  ranges  using  various  conditions.  Option  6  prints  values  of  basic  wind 
statistics  using  a  rotated  coordinate  system,  and  Option  8  prints  wind  speed  percentiles.  Option  #9 
generates  a  table  of  probabilities  of  ranges  of  wind  directions  and  speeds. 

3.3.1  Output  Formats.  The  format  of  output  varies  according  to  the  option  being  used.  The 
different  formats  are  listed  below. 

a.  Option  #0  prints 
STOP 

END  OF  EXECUTION 
CPU  TIME:  bibb.bb 
EIAPSED  TIME:  e:ee.ee 

where  b:bb.bb  is  the  CPU  time  used  during  the  execution  of  the  program,  and  e:ee.ee  is  the  actual 
time  elapsed  since  the  program  was  started.  Both  times  are  in  minutes,  seconds,  and  hundredths  of 
seconds.  A  program  stop  due  to  the  use  of  an  illegal  option  number  prints  the  same  type  of  output. 

o.  The  output  produced  by  option  #1  begins  with  a  header,  "X  YLOW  YHIGH"  and  then  continues 
with  lines  formatted  in  such  a  way  that  the  appropriate  information  is  positioned  below  each  header 
word.  The  output  lines  have  the  format,  XXXX.XXXX.YYYY.YYYY  ZZZZ.ZZZZ  where  the  values  XX...., 
YY....,  and  ZZ....  correspond  to  xl,  yl,  and  y2  where  (xl,yl)  and  (xl,y2)  are  the  intersection 
points  of  the  line  x  =  xl  with  the  ellipse  to  be  plotted.  This  setup  is  shewn  graphically  in  Figure 
3. 


Figure  3.  Plotting  Points  Output  by  Option  #1. 
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One  set  of  points  is  printed  per  line,  so  that  the  number  of  lines  printed  with  this  fonrat  i6  equal 
to  the  number  of  point  sets  specified  by  the  user.  Ihe  reference  information  printed  at  the  er*d  of 
this  option's  output  consists  of: 

EXTREME  VALUES  OF  X:  xxxx.xxxx  xxxx.xxxx 
EXTREME  VALUES  OF  Y:  yyyy.yyyy  yyyy .yyyy 

where  xxxx.xxxx  and  yyyy.yyyy  are  the  largest  and  smallest  values  of  X  and  Y  on  the  elliptical 
contour . 


c.  Option  12  prints  a  single  line  of  output  with  the  format 
PROBABILITY  OF  GIVEN  WIM)  DIRECTION  RANGE  .pppp 

where  .pppp  is  the  probability  asked  for  by  the  user  (always  between  0  and  1). 

d.  Option  #3  prints  a  single  line  of  output  with  the  format 
PROBABILITY  OF  GIVEN  RANGE  OF  SPEEDS  .pppp 

where  .pppp  is  the  probability  of  the  range  of  speeds  given. 

e.  Option  #4  prints  a  single  line  of  output  with  the  format 
PROBABILITY  OF  GIVEN  RANGE  OF  SPEEDS  AND  DIRECTNS.  .pppp 

where  .pppp  is  the  joint  probability  of  the  range  of  directions  and  speeds  given  by  the  user. 

f.  Option  #5  prints  a  single  line  of  output  with  the  format 
PROBABILITY  OF  GIVEN  RANGE  OF  SPEEDS  AT  GIVEN  DIRECT  .pppp 

where  .pppp  is  the  probability  of  the  speed  range  specified  when  the  wind  direction  is  equal  to  the 
specified  value. 

g.  Option  #6  prints  a  single  line  of  output  with  the  format 
NEW  STDEVX:  xxxx.xx  NEW  STDEVY:  yyyy.yy 

where  xxxx.xx  and  yyyy.yy  are  the  standard  deviations  of  the  zonal  and  meridional  coordinates  in  the 
rotated  coordinate  system.  The  option  prints 

NEW  MEAN  X:  iiiimn . rrm  NEW  MEAN  Y:  nnnn.nn  NEW  CORR:  C.CCCC 

where  mrrmm.mm  and  nnnn.nn  are  the  means  of  the  zonal  and  meridional  wind  components  in  the  rotated 
coordinate  system,  and  c.cccc  is  the  new  correlation  coefficient. 

h.  Cation  #7  does  not  print  any  data. 

i.  Option  #8  prints  a  single  line  of  output  with  the  format 
CORRESPONDING  WIND  SPEED:  pppp.pp 

where  pppp.pp  is  the  percentile  value  of  the  wind  speed  asked  for  by  the  user. 

j.  Option  #9  prints  one  page  of  output,  containing  a  table  of  probabilities  of  various  ranges 
of  wind  directions  and  speeds.  An  example  is  shown  in  Appendix  A. 

3.3.2  Sample  Output.  Appendix  A  shows  an  example  of  program  usage  in  which  an  analyst  uses  every 
possible  option. "Details  on  this  example  are  given  in  that  section  of  the  manual. 

3.4  Utilization  of  System  Outputs.  The  output  produced  by  this  program  is  typically  used  to  answer 
customer  requests  for  probabilities  of  wind  direction  intervals  and/or  wind  speed  ranges  or  wind 
speed  percentiles. 
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j.*»  Recovery  and  Error  Correct  ion  Procedures 

No  explicit  restart  capability  is  provided.  The  standard  BBC  detaching  capabilities  apply:  a  user 
can  sometimes  continue  running  the  prograei  from  where  it  was  interrupted  if  the  interruption  woe  due 
to  a  connect  ion  failure,  local  system  problems,  etc.  The  user  can  also  explicitly  abort  the  program 
using  a  control  *Ca  if  repeated  errors  occur  or  if  the  program  is  running  extremely  slowly  due  to 
heavy  machine  usage,  and  then  start  again  frcm  the  beginning. 

a  user  can  ’detach*  the  program  if  he  is  using  option  #9  and  sending  his  output  to  a  data  file. 
This  procsAire  is  recommended  in  cases  where  machine  usage  is  substantial  and  this  option  My 
require  up  to  a  half  hour  to  complete  its  output.  The  user  should  than  reattach  to  the  program  at  a 
later  time,  preferably  about  a  half  hour  after  detaching  from  it.  Option  *9  will  then  indicate  if 
its  output  to  the  data  file  is  complete.  Detaching  the  program  while  using  any  other  option  causes 
the  program  to  suspend  execution  and  wait  until  the  user  reattaches  to  it. 

Erroneous  basic  wind  statistics  can  be  corrected  by  using  option  #7  to  input  the  correct  set  of 
values.  The  original  basic  statistics  can  also  be  retrieved  if  a  faulty  change  of  coordinates 
occurs  in  option  #6.  This  is  done  by  using  that  option  again  with  an  angle  of  rotation  equal  to  the 
previous  value  input  multiplied  by  -1.  The  user  can  run  any  other  option  repeatedly  using  different 
input  values  if  its  initial  results  were  wrong  because  of  erroneous  input. 


Appendix  A 

SAMPLE  BUN  OF  TOE  HIND  STATISTICS  PKX3UM 


The  ranters  in  the  following  description  correspond  to  the  sections  indicated  on  the  sample  run. 

1.  The  user  begins  by  starting  the  program  with  DNMST.  The  program  then  prints  at  the  user's 
request,  a  list  of  the  available  output  options.  It  then  asks  the  user  to  input  a  set  of  basic  wind 
statistics. 

2.  The  user  decides  to  get  IS  points  on  the  55«  probability  contour  and  invokes  option  II. 
Note  that  the  first  and  last  point  set  are  on  the  left  and  ri^it  tips  of  the  ellipse;  their  values 
of  YIjQW  and  THIGH  are  identical. 

3.  Option  12  is  used  to  get  the  probability  of  a  wind  direction  range  from  270°  (W)  through  N 
to  30°  (WE) . 

4.  Option  13  is  used  to  get  the  probability  of  a  wind  speed  range  from  10  to  30  wind  speed 
units. 

5.  The  user  selects  option  14  and  gets  the  joint  probability  of  a  wind  speed  between  10  and  2S 
units  and  a  wind  direction  between  100  and  200  degrees. 

6.  Option  IS  gives  the  user  the  probability  of  the  wind  speed  being  between  5  and  SO  units, 
assuming  a  south  wind  (190°). 

7.  It's  time  for  a  20  degree  counterclockwise  rotation  of  coordinates.  Option  #6  does  this  and 
prints  new  basic  statistics. 

8.  The  user  invokes  option  17  and  inputs  a  new  set  of  basic  wind  statistics. 

9.  Option  18  is  used  to  get  the  88th  wind  speed  percentile.  The  wind  speed  will  be  below  this 
value  88  percent  of  the  time. 

10.  The  user  generates  a  table  of  probabilities  for  various  ranges  of  wind  directions  and 
speeds  using  option  19.  The  location  involved  was  Andrews  Are  and  the  period  of  record  used  was 
January  1960  to  Decanter  1974.  The  wind  speed  units  used  were  miles  per  hour.  Note  that  the  total 
probability  figure  is  not  exactly  1  (in  this  case,  the  result  is  1.0001)  due  to  truncation  and  round 
off  errors. 

11.  The  user  steps  the  program  by  using  option  10.  Total  CPU  tune  used  was  about  3  minutes  and 
47  seconds  and  the  session  lasted  about  17  minutes. 
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Appendix  H 

WIND  srAHSIlCS  MODI  I  FLOW  CHART 
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PROCRAM  It>  -  1)N  DWKT 

programmer  and  analyst  -  2LT  novocrad.  usafetac/dnd 

CREATED  ON  -  15  MAY  79  FOR  PROJECT  1893 
VERSION  -  1.0 

DESCRIPTION  -  THIS  PROGRAM  INTERACTIVELY  CALCULATES  A  VARIETY 
OF  WIND  RELATED  STATISTICS  AND  PROBABILITIES,  USING  AN  ASSUMED 
BIVARIATE  NORMAL  DISTRIBUTION  OF  THE  ZONAL  AND  MERIDIONAL 
KIND  COMPONENTS .  ITS  DATA  BASE  CONS  ISIS  OF  A  SET  OF  FIVE 
BASIC  STATISTICS  WHICH  DESCRIBE  HIE  WIND  AT  A  CIVEN  ALTITUDE, 
LOCATION  AND  CLIMATOLOGICAL  TIME  PERIOD.  THESE  BASIC  STATISTICS 
ARE  THE  MEANS  AND  STANDARD  DEVIATIONS  OF  THE  ZONAL  AND  MERIDIONAL 
WIND  COMPONENTS ,  AND  THE  1 R  CORRELATION  COEFFICIENT.  THE  USER 
SELECTS  THE  OPTION  WHICH  GENERATES  THE  PARTICULAR  TYPE  OF  WIND 
STATISTICS  DESIRED  BY  TYPING  IN  HIE  APPROPRIATE  OPTION  NUMBER. 

THE  PERMISSIBLE  OPTION  NUMBERS.  AND  THE  TYPES  OF  OUTPUT  WHICH  THEY 
GENERATE,  ARE  AS  FOLLOWS 

OPTION  NUMBER  TYPE  OF  OUTPUT  GENERATED 


1  COORDINATES  OF  POINTS  ON  AN  ELLIPTICAL 
PROBABILITY  CONTOUR 

2  PROBABILITY  OF  A  RANGE  OF 
WIND  DIRECTIONS 

3  PROBABILITY  OF  A  RANGE  OF  WIND  SPEEDS 

4  JOINT  PROBABILITY  OF  A  RANGE  OF  WIND 
SPEEDS  AND  DIRECTIONS 

5  CONDITIONAL  PROBABILITY  OF  A  RANCE  OF 
WIND  SPEEDS  CIVEN  A  WIND  DIRECTION 

6  NEW  BASIC  STATISTICS  USING  A  ROTATED 
COORD  I  NATE  SYSTEM 

7  NO  OUTPUT  GENERATED  -  ALLOWS  USER  TO 
INPUT  NEW  BASIC  STATISTICS  DIRECTLY 

8  PERCENTILE  WIND  SPEEDS 

9  A  TABLE  OF  JOINT  PROBABILITIES 

OF  RANGES  OF  WIND  DIRECTIONS  AND  SPEEDS. 

«  NO  OTITUT  GENERATED  -  STOPS  THF.  PROGRAM. 

THE  PROCRAM  IS  DE8IGNED  TO  RUN  IN  AN  INTERACTIVE  ENVIRONMENT 
AND  ALWAYS  TYPES  OUT  A  REQUEST  FOR  INPUT  WHEN  IT  NEEDS  IT. 

METHOD  -  THE  EQUATIONS  AND  THEOHY  USED  TO  DEVELOP  THIS  PROGRAM 
WERE  TAKEN  FROM  NASA  I'M  X-73319,  "VECTOR  WIND  AND  VECTOR  WIND 
SHEAR  MODELS.  .  .  ",  BY  O.E.  SMITH.  ALL  THE  PAGE 

AND  EQUATION  NUMBERS  GIVEN  IN  THIS  PROGRAM'S  DOCUMENTATION 
REFER  TO  THIS  TECHNICAL  MEMORANDUM.  THE  EQUATIONS  USED  TO  GENERATE 
SETS  OF  POINTS  ON  AN  ELLIPTICAL  PROBABILITY  CONTOUR  (OPTION  *1) 

ARE  LISTED  ON  PP.11-12.  THE  PROBABILITY  OF  A  RANGE  OF  WIND 
DIRECTIONS  IS  COMPUTED  BY  NUMERICAL  INTEGRATION  OF  EQUATION 
(33).  THE  CONDITIONAL  PROBABILITY  OF  A  RANGE  OF  WIND  SPEEDS 
GIVEN  A  WIND  DIRECTION  (OPTION  -5!  IS  CALCULATED  USING  EQUATION 
(30).  THF,  NUMERICAL  INTEGRATION  OF  (36)  OVER  ALL  DIRECTIONS 
IS  USED  TO  COMPUTE  THE  PROBABILITY  OK  A  RANGE  OF  SPEEDS  (OPTION 
#3),  AND  SUCH  AN  INTEGRATION  OVER  A  SPECIFIED  RANGE  OF  DIRECTIONS 
GIVES  THE  JOINT  PROBABILITY  OF  A  RANCE  OF  SPEEDS  AND  DIRECTIONS 
(OPTION  *4).  THE  CALCULATION  OF  NEW  BASIC  STATISTICS  IN  A  ROTATED 
COORDINATE  SYSTEM  (OPTION  *6)  IS  ACCOMPLISHED  USING  EQUATIONS  (45) 
THROUCH  (49).  PERCENTILE  WINDSPEEDS  (OPTION  • 8 )  ARE  COMPUTED  USING 
THE  METHOD  OF  HALF  INTERVALS.  THE  PROBABILITIES  WRITTEN  OUT  BY 
OPTION  #9  ARE  GENERATED  USING  REPEATED  CALLS  TO  THE  SUBROUTINE 
USED  FOR  OPTION  »4. 
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C*  REFERENCES  -  I.  NOVOCRAD,  BENJAMIN  B.  ,  1980  USAFETAC  WIND 

C*  STATISTICS  PIODEI.  I'SEILS  MANUAL,  USAFETAC  TN  80-0«:i. 

C*  2.  SMITH.  O.E..  1970  VECTOR  WIND  AND  VECTOR  WIND  SHEAR 

C*  MODELS  AT  0  TO  27  KM  ALTITUDE  FOR  CAPE  KENNEDY,  FLORIDA. 

C*  AND  VAN DENBERF'G  AFB.  CALIFORNIA. 

C*  1NPIT  -  A  SET  OF  BASIC  WIND  STATISTICS.  PLUS  AN  OPTION  NUMBER  AND  ANY 

C*  ADDITIONAL  INFORMATION  REQUIRED  BY  THE  PARTICULAR  OPTION! 8)  USED.  DETAILS 

C*  ARE  GIVEN  IN  SECTION  0.2  OF  THE  USERS  MANUAL. 

C*  OUTPUT  -  THE  STATISTICS  CALCULATED  BY  THE  PARTICULAR  OPTION  SELECTED  BY 

C*  THE  USER.  DETAILS  ARE  GIVEN  IN  SECTION  3.3  OF  THE  USERS  MANUAL. 

C*  SYSTEM  SUBPKOGAMS  USED  -  COS.  SIN.  SORT 
C** 

C*  USER  SUBPROGRAMS  USED  -  A.  B,  C.  D.  EI.IPSE,  PDIR.  PNORM,  PRDIR,  PR8RD, 

C*  ROTATE,  RSPGDR,  SLEGDR,  SPCNT,  TRANS.  WROSE. 

C*  EST I  MATED  CPU  TIME  -  THE  CPU  TIME  USED  IS  VARIABLE  AND  DEPENDS  ON 
C*  WHICH  OUTPUT  OPTION! S)  ARE  SPECIFIED  BY  THE  USER.  DETAILS  ARE 

C*  GIVF.N  IN  SECTION  2.5.3.  OF  THE  USERS  MANUAL. 

C** 

C************ ******************************************************** 

ccc 

c  ***  MAIN  PROGRAM  *** 

CCC 

DIMENSION  ILOCNI5) ,  1 LVLYR!  6 ) , IPOR<4) . IWSUI2) 

COMMON  /Bl.K!/  SQSTX.SQSTY,  PROSIT) 

COMMON  /BLK2/  XBAR.YBAR 
COMMON  /BLK3/  CORH.  UF.NOM 
DATA  I TERM,  I  YES/  - 1 . ’ Y  '/ 

I  FORMAT  !*  ***  USAFETAC/DND  WIND  STATISTICS  PROGRAM  ***’,/ 

CCC 

C  **  ASK  USER  IF  INFORMATION  ABOUT  OUTPUT  OPTIONS  IS  NEEDED  ** 

CCC 

S’  DO  YOU  NEED  INFORMATION  ABOUT  VARIOUS  OUTPUT  OPTIONS? 

•  •  TYPE  A  ”  Y”  IF  YES,  A  ’’ll”  IF  NO’) 

WRITE  !  I TERM ,  I) 

62  FORMAT!  A4) 

READ  !  I TERM, 62)  I  HELP 

IF  !  IHELP.NE.  IYES)  CO  IX)  70 

CCC 

c  ***  THIS  SECTION  HELPS  OUT  A  NEW  USER  *** 

CCC 

WRITE  !  I TERM , 4 ) 

WRITE  !  I TERM, 44) 

4  FORMAT  !’  ***  OUTPUT  OPTIONS  AVAILABLE  ***’,/ 

S’  0  STOPS  THE  PROCHAM’,/ 

S’  1  PRINTS  POINTS  DEFINING  AN  ELLIPTICAL  PROBABILITY  CONTOUR’,/ 

S’  2  PRINTS  OUT  PROBABILITY  OF  THE  WIND  DIRECTION  BEINC  WITHIN’,/ 

S’  A  SPECIFIED  RANGE’,/ 

S’  3  PRINTS  OUT  PROBABILITY  OF  THE  WIND  SPEED  BEING  WITHIN  OR’,/ 

S’  OUTSIDE  OF  A  SPECIFIED  RANGE’ , / 

S’  4  PRINT’S  OUT  PROBABILITY  THAT  THE  WIND  SPEED  IS  WITHIN’,/ 

S’  OR  OUTS  I DE  OF  A  SPECIFIED  RANGE  AND  THAT  THE  WIND  DIRECTION’,/ 

S’  IS  WITHIN  A  SPECIFIED  RANGE’) 

44  FORMAT  !’  5  PRINT'S  OUT  THE  (CONDITIONAL)  PROBABILITY  OF’,/ 

S’  THE  WIND  SPEED  BEINC  WITHIN  OR  OUTSIDE  OF  A  SPECIFIED’./ 

S’  RANGE  WHEN  THE  WIND  DIRECTION  IS  KNOWN’,/ 

S’  6  ROTATES  THE  COORDINATE  AXES  THROUGH  A  SPECIFIED  ANGLE  AND  ’./ 

S’  POINTS  OUT  THE  RESULTING  VALUES  FOR  THE  5  BASIC  PARAMETERS’,/ 

S’  7  INPUTS  NEW  BASIC  PARAMETERS’,/ 

S’  8  PRINTS  OUT  A  PERCENTILE  BOUNDARY  VALUE  FOR  THE  WINDSPEED’,/ 

S’  9  PRINTS  OUT  A  PROBABILITY  TABLE  FOR  THE  WIND  DIRECTION  AND’, 

S’  SPEED’) 

CCC 

c  ***  INPUT  THE  FIVE  BASIC  PARAMETERS  *** 

CCC 

70  WRITE  !  I TERM, 77) 

77  FORMAT  !’  INPUT  MEAN  X , STDE VX , ME AN  Y, STDEVY, CORK.  COEFF.’> 

READ  ! I TERM. 88)  XBAR, STDEVX, YBAR. STDEVY, CORR 
88  FORMAT  !  RF) 

IF  (STDEVX. GT.0. AND. STDEVY. GT.0. )  GO  TO  301 
WRITE  (  I TERM , 300 ) 

300  FORMAT  (  ’  STDEVX  AND  STDEVY  MUST  BE  GREATER  THAN  ZERO’ ./ 

S’  PLEASE  TYPE  IN  THE  BASIC  PARAMETERS  AGAIN’) 

CO  TO  70 

301  IF  (  CORF. CT. -  I . AND. CORR. LT.  1 . )  GO  TO  303 
WRITE  (  I  TERM, 302) 

302  FORMAT  (’  THE  CORRELATION  COEFFICIENT  MUST  BE  BETWEEN  -1  AND  I.',/ 

S’  PLEASE  TYPE  IN  THE  BASIC  PARAMETERS  AGAIN.’) 

GO  TO  70 
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cco 

C  **  -  CALCULATE  RECURRING  TERMS  SAVED  IN  COMMON  STORAGE  *»* 

CCC 

Don  SOSTY  “  STDEVY  *  STDKVY 
SOSTX  --  STDEVX  *  STDEVX 
DF.NOM  =-  I  -  (CORR»GORR) 

PROSTD  -  STDEVX  *  STDEVY 
ccc 

C  ***  ASK  DSF.il  WHAT  OPTION  HE  WANTS  *** 

CCC 

5  WHITE  (  ITERM. 95 > 

95  KORMATl  ’  INPUT  DESIRED  OPTION’) 

READ  (  ITERM.97)  lOPTN 
97  FORMAT  (115) 

CCC 

c  ***  WARN  THE  USER  IF  A  NONEXISTENT  OPTION  IS  USED  *** 

CCC 

IF  (  IOPTN.  IT  .0. OR. IOPTN. GT. 9)  GO  TO  990 

CCC 

c  ***  GO  TO  THE  SECTION  GOVERNING  THE  APPROPRIATE  OPTION  *** 

CCC 

IF  ( IOPTN. EQ.0)  GO  TO  999 

GO  TO  ( 10,20,30,40,50.60,70.80,90) , IOPTN 

CCC 

c'  ***  THIS  SECTIONS  CONTROLS  OPTION  *1  *** 

CCC 

10  WRITE  (  I TERM,  1 10) 

110  FORMAT  (’  INPUT  PROBABILITY  VALUE’) 

READ  ( ITERM, 1 12)  PROD 

112  FORMAT  (IF) 

IF  (PROB.LE.0. .OR. PROD. GE. 1.)  GO  TO  990 

113  FORMAT  (  ’  INPUT  NUMBER  OF  X-YLOW-YHICH  SETS  DESIRED’) 

WRITE  (  ITERM,  1 13) 

READ  (ITERM, 97)  NPTS 
WRITE  (  (TERM.  Ill) 

111  FORMAT  (  5X.  ’X’ , 8X,  ’ YLOW’ , 6X,  ’ YHIGH’ . /3(  IX, 9  (  *-’))) 

CCC 

c  ***  USE  “EL IPSE"  TO  CALCULATE  AND  WRITE  OUT  REQUIRED  POINT  SETS  *** 
CCC 

CALL  EL  IPSE  (  PROB, X, Y, XSMAIX. XLARCE. Y8MALL, YLARGE, NPTS) 

CCC 

c  ***  PRINT  OUT  EXTREME  VALUES  OF  X  AND  Y  *** 

CCC 

WRITE  (  ITERM, 1 14)  XS MALL, XLARCE 
WRITE  (ITERM.  116)  YS  MALI.,  YLARGE 
I  14  FORMAT  (IX,’  EXTREME  VALUES  Of  X  ’ , 2( IX, F9. 4> ) 

116  FORMAT  (IX.’  EXTREME  VALUES  OF  Y  ’ , 2(  1 X, F9 . 4) ) 

GO  TO  5 

CCC 

C  ***  THIS  SECTION  CONTROLS  OPTION  *2  *** 

CCC 

20  WRITE  (  ITERM,  120) 

120  FORMAT  (’  INPUT  RANGE  OF  DIRECTIONS’) 

READ  (  ITERM.  122)  DIRI.DIR2 

IF  (DIR1.NE.DIR2)  CO  TO  21 
WRITE  (  ITERM. 121) 

121  FORMAT  (  ’  THE  TWO  INPUT  DIRECTIONS  CAN’ ’T  BE  THE  SAME.’) 

CO  TO  20 

21  DIRl  =  TRANS(DIRI) 

DIR2  =  TRANS(DIR2) 

122  FORMAT  (2F) 

P  =  PRDIR(UIRI,DIR2) 

IF  (P.LT. l.E-04)  P  =  0. 

WRITE  (  ITERM.  124)  P 

124  FORMAT  (’  PROBABILITY  OF  GIVEN  WIND  DIRECTION  RANGE’, F6. 4) 

GO  TO  5 

CCC 

C  ***  THIS  SECTION  CONTROLS  OPTION  *3  *** 

CCC 

30  WRITE  (  ITERM,  130) 

130  FORMAT  (’  INPUT  RANGE  OF  SPEEDS  ’) 

READ  ( ITERM, 132)  SPDI , SPD2 

132  FORMAT  (2F> 

IF  (SPDI. GE.0. AND. SPD2.GE.0.)  GO  TO  32 
WRITE  (  ITERM,  131) 

131  FORMAT  C  ONE  OR  MORE  OF  THE  INPUT  SPEEDS  ARE  NEGATIVE* ,  / 

*’  PLEASE  INPUT  THE  SPEEDS  AGAIN.’) 

GO  TO  30 

32  CALL  PRSIUX0.  ,360.  .SPDI,  SPD2.  PROB) 

IF  (  PROD. LT.  1 . E-04)  PROB  =  0. 

WRITE  (  ITERM, 134)  PROB 

134  FORMAT  (’  PROBABILITY  OF  GIVEN  RANGE  OF  SPEEDS' , F6 . 4) 

CO  TO  5 
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coc 

C  ***  THIS  SECTION  GONTHOLS  OrTlON  *4  *** 

COC 

40  WHITE  (  I TERM, 120) 

HEAD  (  ITERM,  122*  l>!KI,DIR2 
IK  (  Dim  .  NE.  1)1  H2>  CO  TO  41 
WRITE  <  1TERM, 121 ) 

CO  TO  40 

41  l)  1  R 1  •'  TRANS)  I)  1 R I ) 

DIR2  =  TRANS) DIR2) 

WRITE  (  I TERM ,  I30> 

READ  (  I  TERM,  132>  Sl’D  1  ,  SPD2 

IK  ( SPD 1 . CE . O . ANP . SPD2 . CE . 0 . )  CO  TO  42 

write:  (  [term,  i:t  1 1 

CO  TO  40 

42  IF  (DIR2.CE.DIR1)  CO  TO  49 

CALL  PRSRD  ( D I R  1 , 360 . , SPD I , 8PD2 , PROB 1 ) 

CALL  PRSRD  (  O  ,  I)  I R2  ,  SPD  1  ,  SPD2 ,  PROB2) 

PROB  =  PROB1  -*•  PR0B2 
CO  TO  5 1 

49  CALL  PRSRD  ( D I Rl . DI R2 . SPD I . SPD2 , PROB) 

51  IF  (PROB.LT. I . F-04)  PROB  =  0. 

WH1TF.  (  1TF.RM.  142)  PROB 

142  FORMAT  ('  PROBABILITY  OF  GIVEN  RANGE  OF  SPEEDS  AMD  DIRECTNS’ ,F7.4) 
GO  TO  5 

ccc 

c  ***  THIS  SECTION  CONTROLS  OPTION  *3  *** 

ccc 

50  WRITE  (  ITF.RM,  150) 

150  FORMAT  ('  INPUT  DIRECTION  ') 

READ  ( ITERM, 112)  DIR 
DIR  =  TRANS) DIR) 

WRITE  !  ITERM,  130) 

READ  ( ITERM, 132)  SPD  1 , 8PD2 

IF  (SPl)l  .  GE.0.  AND.SPD2.GE.  0.  >  CO  TO  52 

WRITE  (  ITERM,  13!) 

GO  TO  50 

32  CALL  RSPGDR!  DIR, SPD 1 , SPD2 , PROB) 

WRITE  (  ITERM,  152)  PROB 

152  FORMAT  (•  PROBABILITY  OF  GIVEN  SPEED  RANGE  AT  CIVEN  DIRECT. F6 . 4) 
CO  TO  5 

CCC 

C  ***  THIS  SECTION  CONTROLS  OPTION  *6  *** 

CCC 

60  WRITE  (  ITERM,  160) 

160  FORMAT  (’  INPUT  ANGLE  OF  ROTATION') 

READ  (ITERM, 112)  ALPHA 

CALL  ROTATE) ALPHA, STDEVX, STDEVY) 

WRITE  I  ITERM,  162)  STDEVX . STDEVY 

162  FORMAT  (’  NEW  STDEVX  ’,F7.2.'  NEW  STDEVY  ’,F7.2> 

WRITE  ( ITERM, 163)  XDAR, YBAR.CORR 

163  FORMAT  ('  NEW  MEAN  X  ’.F7.2,’  NEW  MEAN  Y  ' , F7.2, '  NEW  CORR  \F6.4) 
GO  TO  5 

CCC 

c  ***  THIS  SECTION  CONTROLS  OPTION  *8  *** 

CCC 

80  WRITE  (  ITERM. 81) 

81  FORMAT  (’  INrUT  THE  REQUIRED  PERCENTILE') 

84  READ  (ITERM, 112)  PCNT 

IF  (PCNT. GT.0. AND. PCNT. LT. 100. >  GO  TO  B6 
WRITE  (ITERM, 89) 

89  FORMAT  (  ’  PERCENTILES  MUST  BE  BETWEEN  0  AND  100.’,^ 

*’  PLEASE  INPUT  THE  NUMBER  AGAIN.’) 

GO  TO  84 

86  PCNT  =  PCNT  /  100. 

CALL  8 PCNT  (  PCNT, PSPD) 

WRITE  (ITERM, 82)  PSPD 

82  FORMAT  (’  CORRESPONDING  WIND  SPEED  ’,F8.4> 

GO  TO  5 

CCC 

C  ***  THIS  SECTION  CONTROLS  OPTION  #9  *** 
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ccc 


90 

WRITE  i 

ITERM, 190) 

190 

FORMAT 

’  INPUT  THE.  PERIOD  OF  RECORD 

-  UP  TO 

20  CHARACTERS’ ) 

READ  «  1 

HM.  191  •  IPOR 

19  1 

FORMAT 

‘  4A5  ) 

92 

WRITE  ( 

1  PERM.  192) 

192 

FORMAT 

(  '  INPUT  HIE  LEVEL  OR  LAYER  - 

UP  TO  30 

CHARACTERS ’ ) 

READ  t  1 

TERM, 193)  ILVLYR 

193 

FORMAT 

(  ft  AS) 

94 

WHITE  i 

ITERM. 1 94 i 

194 

FORMA E 

■  INPUT  T1IF.  LOCATION  -  UP  TO 

25  CHARACTERS  ’  ) 

READ  t  1 

PERM.  D5>  II.OCN 

195 

FORMAT' 

<  :>  vr. ' 

96 

WRITE  ( 

ITERM. 196) 

196 

FORMAT 

t  •  INPUT  THE  WIND  SPEED  UNITS 

-  UP  TO 

10  CHARACTERS’ ) 

READ  (  ITERM.  197)  I WSU 

197 

FORMAT 

l  1!  ) 

CALI.  WHOSE.  1  IPOR,  ILOCN.  ILVLYR,  IWSU) 

GO  TO  5 

CCC 

c  ***  PROGRAM  ABORT  FOR  A  NONEXISTENT  OPTION 

*** 

CCC 

990  WRITE  t  i TERM. 900) 

900  FORMAT  I  ’  I  I.LF.GAJ,  OPTION  NUMBER  GIVEN  -  PROGRAM  STOPPED.’) 
999  STOP 
END 

CCG 

c  ***  FUNCTION  A  -  COMPUTES  TERM  "A”  P.23  *** 

ccc 

FUNCTION  A( S I ND I H, COSD 1 R, CORK, DENOM, N) 

COMMON  /BLKl/  SQSTX, SQSTY, PROSTD 

A  =  (  (  COSI)in*COSDIR/SOSTY)  -  ( 2*C0RR*C0SDIR*S INDI R/PROSTD) 
*+  (  S  I  ND  I  R*S  INDI  R/SQSTX)  )  /  DENOM 
IF  (N.EQ.2)  RETURN 
A  =■  SORT!  A) 

RETURN 

END 

CCC 

c'J  ***  FUNCTION  R  -  COMPUTES  TERM  “B"  P.23  *** 

CCC 

FUNCTION  R(  S  I  Nl)  I R,  COSD  I R,  CORR,  DENOM) 

COMMON  /ISLK1/  SQSTX, SQSTY, PROSTD 
COMMON  /BLK2/  XBAR.YBAR 

H  =  (  < YBARXCOSDUI/SQSTY)  -  ( CORR* ( YBAR*S I HD I R+  XBAR*COSD I R) 
*  /PROSIT)'  +  ( XBAR*S INDI R/SQSTX) > /DENOM 
RETURN 
END 

CCC 

C  ***  FUNCTION  C  -  COMPUTES  TERM  "C”  P.23  *** 

CCC 

FUNCTION  C( CORR, DENOM. N) 

COMMON  /BLKl  /  SQSTX,  SQSTY,  PROSTD 
COMMON  /BL.K2/  XBAR.YBAR 

C  =  ( < YBAR*YBAR/SQSTY)  -  ( 2*CORR*XBAR*YBAR/PROSTD) 

*+  ( XBAR/XBAR/SQSTX) ) /DENOM 
IF  (N.EO.2)  RETURN 
C  =  SORT'D 
RETURN 
END 

CCC 

C  ***  FUNCTION  I)  -  COMPUTES  TERM  "D"  P .  23  *** 

CCC 

FUNCTION  D( DENOM) 

COMMON  /Bl.Kl/  SQSTX,  SQSTY,  PROSTD 
DATA  TWOF 1/6. 2831 852/ 

1)  =  1  /  (  TWOP  I  /-PROSTD/ SORT!  DENOM)  ) 

RETURN 

END 

CCC 

C  ***  SUBROUTINE  EL  IPSE  -  CALCULATES  AND  PRINTS  POINTS  ON  AN 

C  ***  ELLIPTICAL  PROBABILITY  CONTOUR 
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m 


*** 

*** 


ccc 


SUBROUTINE  i  t  I  rst'.t  PROB.X,  V,  X8HAL1..  XLAJICE,  YSMALL,  YLAHCE,  NPT9) 

COMMON  BLK1  SOSTX . SQSTY ,  PROSTR 

COMMON  -i J.iCT  M<AR.  YDAR 

COMMON  'BUCll  CoHR, DENOM 

REAL  LMBD.  l.MRDE,  LYROOT 

o  -•  ai.oo  i  -  prob> 

LMBD  =  SORT (  -2 .  *  DENOM*Q) 
l.MBOF  -  SORT'  2.  *0) 

A  -  SOSTY 
IT  =  -  I 

B  =  -2.  *CORR*TROSTD 

c  -  sostx 

0  -  -  1  .  *(  (  B*VHAR)+<  2.*A*XHAR>  ) 

E  =  I  s- •  «  B* XB  All)  +  (  2  .  *C*YBAR)  ) 

F  =  t  ,\*  Xh  Ui*XUMU  t  (  C*YB\R*YBAl\)  +  <  B*XBAR*YBAR> 

*  -  (  .A*C*1.MBD*LMBD' 

XV  ACT  -■  LfIBOK  *  SORT!  SOSTX) 

YKACT  l.MBBF  *  SORTl SQSTY) 

XSN.ALi.  =  XBAU  -  XT  ACT 

XLARCF.  =  XB,\H  -r  XFACT 

YSMALL  -  YU  AM  -  YF  ACT 

YLAROE  =  YBAR  +  YKACT 

XDIKK  =  XI.ARCE  -  XSMALL 

NPTS  =  N>“  I  S  -  l 

DO  20  N  -  I.NPTSfl 

XTKMP  =  XSMALL  3  XDIFK*(  N-  l ) /NPTS 

QUADA  -  C 

OUADB  -  l  R*XTEMP>  +  E 

Q1IADC  -  ( A*XTEnP*XrEHP>  +  <  D*XTEMP)  +  F 
QROOT  =  <QUADB*QUA1)B)  -  ( 4 .  *QUADA*OUADC) 

IF  (  QROOT.  LT,  O)  QROOT  -  0. 

OF  ACT  =  SORTl  QROOT) 

LYROOT  =  ii-l.  *U!JAI>8)  -  SQRT(  QROOT) ) /(  2.  #QUADA> 

UY1100T  =  (  (  -  1  .*QUADB)  +  SORTC  QROOT)  )  /  (  2 .  *QUADA) 

222  FORMAT  Ci(  IX.F9.4)) 

WHITE  (IT,  222 '  XTEMP ,  LYROOT ,  UYROOT 
20  CONTINUE 
RETURN 
END 

CCC 

C  **  FUNCTION  PD  I II  -  OiVES  PROBABILITY  DENSITY  FOR  WIND  DIRECTION  ** 
CCC 

FUNCTION  PD 1 R( TRETA) 

COMMON  /BLKI/  SQS  TX.  SQSTY,  PROSTD 
COMMON  /-IJLK2.  XBAll.YBAR 
COMMON  /RLK3/  COHH,  DENOM 
DATA  SQRT2IV2. 5001.28/ 

S1NDIR  =  S  INI  THETA) 

COSDIR  =  COS< THETA) 

ASO  =  A( S I NDI R. CONI) I R. CORR, DENOM,  2) 

A I  =  SORT!  ASO' 

B1  *  R(S1NP!R. COSDIR, CORR, DENOM) 

BOV A  -  B I  /  A I 
BOVASO  =  BOV A  *  HOVA 

PDIR  =  (  D(  DENOM) /ASO)  *  EXP(  - . 5#C<  CORR, DENOM, 2) )  * 

•  (  I  +  (  SQHT2P*B0VA*EXP( . 3*B0VASQ) *PNORM( BOVA) ) ) 

RETURN 
END 


FUNCTION  PRDiR  -  COMPUTES  PROBABILITY  OF  WIND  DIRECTION  RANCE  ** 
FUNCTION  PRDIRl  D1R1 ,DIR2) 

DATA  CONV/I . 745329E-02/ , PRDIR, PRD IRO/2*0 .  ITERM/-1/ 

NSLICE  =  20 
PRDIR  =  0. 

IF  <  D1R2  -  DIRI>  10,110,30 
SLSIZF.  =  (  360.  -DIUl ) /NSLICE 
DO  12  IC  =  I .NSLICE 

ZD  I R  =  DIRI  +  ( IC*SLNIZE)  -  (SLSIZE/2.) 

PRDIR  =  PRDIR  +  ( PDIRl ZD1R*C0NV) ) *SL8IZE*C0NV 
SLSJZE  =  DIR2/NSLJCE 
DO  14  IC  =  l, NSLICE 
ZD  1 II  =  (  IC.tSI.SIZF.)  -  (SLSIZE/2.) 

PRDIR  =  PRDIII  +  <  l‘DIR<  ZDIR*C()NV)  >  *SLSlZE*CONV 
IF  (PRDIR  -  PRDIRO.I.T.  l.E-05)  RETURN 
NSLICE  NSLICE  *  2 
PRD  I RO  =  PRDIR 
CO  TO  5 


CCC 
C  ** 
CCC 


5 

10 


12 


14 
I  b 


n 


30  Sl.Si/I.  =  < DIR2-DIR1)/NSLICE 
DO  3-J  1C  =  l.NSLICE 

Zl)  1 H  mm  +  »  IC*SLSIZE>  -  ISL81ZE/2. > 

32  PRD  in  =  1*RD  I  It  +  (  PDIIU  ZOIR*CONV>  > *SLS lZE*CONV 
GO  I'D  If) 

END 

ccc 

i **  SUBROUTINE  PRSK!)  -  COMPUTES  THE  PROBABILITY  OF  A  RANGE  OF  ** 

C  **  SPEEDS  AND  A  RANGE  OF  DIRECTIONS  ** 

CCC 

SUBROUTINE  PHSRDC  D1R1, D1R2.SPD1 .SPD2.PR0B) 

DATA  GONV-  1 . 743U2DK-02/ 

I  TERR  =  -l 
,JSl.  ICE  -  2D. 

5  PROD  =  D. 

IK  i UIH2.E0. D1RI>  CO  TO  20 
SI.SIZE  i  DIR2-DIRD  AJSLICE 
DO  10  IN  =  l.JSI.ICE 

/.Dili  DIR1  +  (  1N*SLSIZE)  -  (SLSIZF./2.) 

CALL  RSPGDH*  ZDIR.SPDI  ,  SPD2 .  PRB I  NT) 

DD2  =  DIHl  +  IN*8LSIZF. 

DDt  DIIU  *  i  IN-  1  >*SLS1ZE 
10  PROB  =  PROB  +  PRB I NT  *  PRD  I R( DD 1 , DD2) 

IF  (  PP.OB  PHOKO.  LT.  1  .  E-04)  RETURN 
.JSI.ICF.  -  JSI.ICK  *  2 
PROBO  =  PROB 
GO  ro  5 

20  CALL  RSPGDRI  Dim.  SPD 1  ,  SPD2 ,  PRB  1  NT) 

PROB  -  PRB  I  NT  *  P1U)1R(  DIR1 ,  D1R1 ) 

return 

END 


C  **  SUBROUTINE  ROTATE  -  ROTATES  COORDINATES  THROUGH  A  CIVEH  ANGLE  ** 
CCC 

SUBROUTINE.  ROTATE!  ALPHA. STDEVX, STDEVY) 

COMMON  /  BLKlz  SQSTX,  SQSTY,  PROSTI) 

COMMON  ZBLK2Z  XBAR. YBAR 
COMMON  ZBLK3/'  COHR.DENOM 
DATA  CON V' 1 . 74332"E-02z 
CVANG  =  ALPHA  *  CONV 
COSA  =  COS! CVANG) 

S1NA  =  SIN!  CVANG) 

COSSQA  -  COSA  *  GOSA 
SINSUA  =  S IN A  *  SIRA 
Ol.DXBH  =  XBAR 
OI.DSQX  =  SQSTX 
OLDSQY  --=  SQSTY 

SQSTX  =  (  SQSTX*COSSQA)  +  ( SQSTY*S I NSQA)  +  C  2*C0RR*PR0STD*C08A*8 1 NA) 
SQSTY  =  (  SQSTY*COSSQA) + ( OLDSQX*S I NSQA)  -  (  2*CORR*PROSTD*COSA*S 1 NA) 
XBAR  =  (XBAR«COSA)  +  ( YBAR*S INA) 

YP.AH  =  < YBAH*COSA)  -  (  OLDXHR*S  I  NA) 

COVXY  =  C.ORR  *  PROSTD 

GOVXY  =  COVXY  *  (COSSQA  -  8INSQA)  +  ( C0SA#8 1 NA*(  OLDSQY-OLDSQX) ) 
STDEVX  =  SORT!  SQSTX) 

STDF.VY  =  SORT!  SQSTY) 

PROSTD  -  STDEVX  *  STDEVY 
CORR  =  GOVXY  /  PROSTD 
RETURN 
END 

CCC 

c  **  SUBROUTINE  RSPGDR  -  COMTUTES  PROBABILITY  OF  A  RANGE  OF  ** 

C  **  STEEDS  GIVEN  A  WIND  DIRECTION  ** 

CCC 

SUBROUT I NE  RSPCDR! D I R, SPD I , SPD2 . PROB) 

IF  (SPD I  .  EQ.0.  )  GO  IT)  10 
IF  (SPD2.EQ.0.)  GO  TO  20 
CALL  SLEGDRI DIR. SPD 1 , PROB 1 ) 

CALL  SLEGDR(  D 1 R, SPD2 , PR0B2) 

PROB  =  PR0B2  -  PROB 1 

IF  (SPD2.1  T.SPD1)  PROB  =  PROB  +  I 

RETURN 

10  CALL  SI.EGDR!  DIR.  SI’D2 ,  PROB) 

30  RETURN 

20  GALL.  SLEGDRI  1)  1  R,  SPD )  ,  PROB) 

PROB  =  I .  -  PROB 

RETURN 

END 
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ccc 

C  **  SUBROUTINE  SLEGDR  -  COMPUTES  THE  PROBABILITY  OF  THE  ** 

C  **  SPEED  BEING  LESS  THAN  OR  EQUAL  TO  A  GIVEN  SPEED  WHEN  ** 

C  **  THE  WIND  DIRECTION  IS  KNOWN  ** 

CCC 

SUBROUTINE  SLEGDR  (  DIR.SPD.PROB) 

COMMON  /BLKl/  SQSTX, SQSTY, PROSTD 

COMMON  /BLK2 '  XBAR, YBAR 

COMMON  /BLK3 '  CORR, DENOM 

DATA  CONV. 1 . 745329E-02/ . SQRT2P/2 . 506628/ 

THETA  -  DIR  *  CONV 
COSTHT  =  COS*  THF.TA) 

SINTHT  =  SIN* THETA) 

ASQ  =  A* SINTHT, COSTHT, CORR, DENOM, 2) 

A 1  =  SORT* ASQ) 

B 1  =  Bt SINTHT, COSTHT, CORR. DENOM) 

BOVA  =  B1  /  A 1 
RS  =  <  Al*SPD)  -  BOVA 
ZX  -  . 5*RS*RS 
IF  (ZX.GT.80. )  ZX  =  80. 

ZZ  =  - . 5*Bl*Bl/ASQ 
IF  (ZZ.LT.-80.)  ZZ  =  -88. 

PROB  =  l  -  (  *  EXP* -ZX)  +  ( SQRT2P*BOVA*< 1 .  -  PNOHM* RS) )) ) 

•  /< EXP*  ZZ)  +  ( SQRT2P*B0VA*PN0RM*  BOVA) ) ) ) 

20  RETURN 
END 


CCC 

C  *  FUNCTION  PNORM  -  GIVES  THE  NORMAL  PROBABILITY  DENSITY  FUNCTION  * 
CCC 

FUNCTION  PNORM*  X) 

DATA  C0/1. 1892071/, Cl/. 23410017/. C2/. 13698952/ , C3/4. 0908724E-4/, 
•C4/2 . 322 1647E-2/ 

IF  <  X. LT.0. )  GO  TO  1 

A  =  ( (  ( C4*X+C3)  *X  +C2)  *X+C1>*X+C0 
IF  (A.GT. I.E20)  A  =  l . E20 
PNORM  =  1 . -  1 . /*  A*A*A*A) 

RETURN 

1  A  =  (  ( ( C4*X-C3)  *X  +C2)  *X  -  Cl)  *  X  +  CO 
IF  (A.GT. 1.E20)  A  =  1 . E20 
PNORM  =  l./(A*A*A*A> 

RETURN 

END 

CCC 

C  **  FUNCTION  TRANS  -  MODIFIES  INPUT  WIND  DIRECTION  SO  THAT  ** 

C  **  THE  WIND  VELOCITY  VECTORS  ALL  8TART  AT  THE  ORIGIN  ** 

CCC 

FUNCTION  TRANS* DIR) 

IF  ( DIR. EQ. 360. )  DIR  =  339.9999 
IF  (DIR.EQ. 180. )  DIR*  179.9999 
IF  (DIR.GE. 180. )  GO  TO  13 
TRANS  *  DIR  +  180. 

RETURN 

15  TRANS  =  DIR  -  180. 

RETURN 

END 

CCC 

c  ***  SUBROUTINE  SPCNT  -  COMPUTES  PERCENTILE  VALUES  OF  VINDSPEED  *** 
CCC 

SUBROUTINE  SPCNT* PCNT, PSPD) 

COMMON  /BLK2/  XBAR. YBAR 
REAL  LOWLIM 
DATA  LOWLIM'0./ 

ULIM  =  SORT*  XBAR* XBAR  +  YBAR* YBAR) /( 1 . -PCNT) 

IF  (PCNT. GT. .83)  ULIM  *  ULIM  /  *-l.*ALOG( l.-PCNT)) 

LOWLIM  =  0. 

5  HLIM  =  (ULIM  +  LOWLIM) /2. 

CALL  PRSRD(0. ,360. ,0. , HLIM, PCNTH) 

IF  (  ABS*  PCNTH-PCNT) . LT. 1 . E-04)  GO  TO  30 
IF  (PCNTH  -  PCNT)  10,30,20 
10  LOWLIM  =  HLIM 
CO  TO  3 

20  ULIM  *  HLIM 
GO  TO  5 

30  PSPD  =  HLIM 
RETURN 
END 
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ccc 

C  **  SUBROUTINE  VROSE  -  CIVES  A  TABLE  OF  FiU)BABILITIES  ** 

C  **  OF  VARIOUS  RANGES  OF  WIND  DIRECTIONS  AND  8PEEDS  ** 

CCC 

SlIBKOUTINL  WHOSE!  IPOR,  ll.OCN,  VLVLYR,  1WSU) 

DIMENSION  1LOCNI 5) .  1LVLYRI  6)  , 1  PORI  4) . RDIRI  16) . RSPIH  I  I) .PROBI  12) 
•  IWSU12) .NMEDIRI  18) .SPDTOTI  12) 

DATA  1TEKM.  1 Y.  1  FORM/  -  1 ,  *  Y  1  ,  "060000000000/ 

DATA  RD11L'  11.25,33.75,56.23.78.75,101.23,123.73,  146.23, 168.73, 
#191 .25.213.75,236.23,258.73,281.25,303.75,326.23.348.73/ 

DATA  HSPD/1 . , 10. ,20. .30. ,40. .50. .60. ,75. , 100. . 130. ,200./ 

DATA  NMED I R/ '  S’,*  SSW  ,  *  SW’  ,  •  WSW‘ ,  *  W  WNW’ , 

#•  NW’  ,  ’  NNV  .  ’  N  ’  ,  ’  NNE’  ,  ’  NF.’  ,  *  ENE’  ,  ’  E  ESE' 

*•  SF.-  ,  ’  SSE’  .  1  CALM1  ,  ’  TOTL'  / 

WRITE  (  1 TERM. 20) 

20  FORMAT  I’  DO  YOU  WANT  THE  PROBABILITY  TABLE  WRITTEN  OUT  TO’, 

#’  A  FILE?  ( Y-  YES , N=  NO) * ) 

READ  l ITERM.23)  1R 

25  FORMAT  I  A5  > 

IF  (  IR. NE.  I Y)  GO  TO  30 
WRITE  (  [TERM. 35) 

35  FORMAT  C ’  PLEASE  INPUT  A  FILENAME  -  UP  TO  3  CHARACTERS') 

READ  (ITERM.23)  IFILE 
WRITE  < I TERM, 26) 

26  FORMAT  (  ’  WOULD  YOU  LIKE  THE  PROGRAM  TO  STOP  AFTERWARDS? ’ ) 

READ  ( I TERM, 25)  IR 

IF  (  IR. EQ.  IY)  ISTOP  =  1 
WRITE  (  ITERM, 27)  IFILE 

27  FORMAT  (  ’  YOUR  PROBABILITY  TABLE  WILL  BE  SENT  OUT  TO’ . 

••  FILE  ’ . A5. ’ .DAT’ ) 

OPEN  (UNIT=20,FILE= IFILE) 

ITERM  =  20 

30  WRITE  (  ITERM. S0>  I FORM, ILOCN, ILVLYR, IPOR, IWSU 

50  FORMAT  <A1./4X,5A5, 1 X. 6A3 , 4A5//. 2X. ’WIND  SPEED  -  ’,2A5/.2X. 

•’DIR  - ’ / , 2X, ’ - */,9X. *1-10  10-20  20-30  30-40  40-80  50-60 

»  60-75  75-100  100-150  130-200  >200  TOTAL’) 

DIRTOT  =  0. 

DO  100  IS  =  1.10 
SPDLOW  =  RSPD< IS) 

SPDHI  =  RSPD( IS+  1 ) 

CALL  PHSRDC 348. 75 , 360 . , SPDLOW, SPDHI , PROBI ) 

CALL  PRSRDi  0 . , 1 1 . 25 , SPDLOW, SPDHI . PR0B2) 

PROB< IS)  =  PROI’,1  +  PROB2 
DIRTOT  =  DIRTOT  +  PROBI IS) 

100  SPDTOTI  IS)  =  PROBI  IS) 

CALL  PRSRD( 348 . 73 , 360 . , 200 . .0. .PROBI) 

CALL  PRSRDI0. . 11.25.200. .0. ,PROB2> 

PROB<  1  I )  =  PROBI  +  PROB2 
PROBI  12)  =  PROBI 11)  +  DIRTOT 
SPDTOTI  12)  =  PROBI  12) 

DO  125  I  =  1. 12 

123  IF  (PROBI 1) -LT. l.E-04)  PROBI I >  =  0. 

WRITE  I  ITERM, 700)  NMEDIRI  1) ,PROB 
700  FORMAT  (  1X,A5.2X,6(F5.4.  IX) ,3<F3.4,2X> ,2(F5.4, IX> ,F6.4/> 

DO  200  ID  =  1  15 
DIRTOT  =  0. 

DO  150  IS  =  1.10 
SPDLOW  =  RSPDI IS) 

SPDHI  =  RSPD< 1S+ I ) 

DIRLOW  -  RDIR( ID) 

DIRH1  =  RDIRI ID+1) 

CALL  PRSRDI DIRLOW, DIRHI, SPDLOW, SPDHI, PROBI  18)) 

DIRTOT  =  DIRTOT  +  PROBI IS) 

130  SPDTOTI IS)  =  SPDTOTI IS)  +  PROBI IS) 

CALL  PRSRDI  DIRLOW,  DIRHI, 200.  .0.  .PROBI  ID) 

PROBI  12)  =  PROBI  11)  +  DIRTOT 
SPDTOTI  12)  =  SPDTOTI 12)  +  PROBI 12) 

DO  175  I  =  1,12 

175  IF  (PROBI  I) .LT. l.E-04)  PROBI  I )  *  O. 

200  WRITE  I ITERM, 700)  NMED I R( ID+1), PROB 
CALL  PRSRDI  0 . . 360 . , 0 . , 1 . , PCALM) 

IF  I PCALM. LT. l.E-04 >  PCALM  =  0. 

WRITE  (ITERM, 710)  NMED IRC! 7) , PCALM 
710  FORMAT  I  IX. A5.2X.F5.4/) 

SPDTOTI 1 )  =  SPDTOTI 1 )  +  PCALM 
SPDTOTI  12)  =  SPDTOTI 12)  +  PCALM 
WRITE!  ITERM, 700)  NMEDIRI  18) .8PDT0T 
IF  I  ITERM. NE. -  1 )  CLOSEIUNIT  *  20) 

ITERM  =  -1 

IF  ( ISTOP. EO. I)  STOP 
WRITE  I  ITERM, 250) 

230  FORMAT  (’  ***  YOUR  PROBABILITY  TABLE  IS  READY  ***’> 

RETURN 

END 
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